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I. INTRODUCTION 


There is a rich literature concerning the transformation properties and the symmetry 
analysis of the general form of the equations that describe the magnetohydrodynamic (MHD) 
equilibrium of a plasma with flows. See e.g., RefiH where transformation properties of a 
generalized Grad-Shafranov (GS) equation are analyzed for cases with and without spatial 
symmetry in a framework that makes it possible to recover both the magnetostatic and the 
“stationary hydromechanic” limits. In Refj2l as an application of the method, solutions 
of the stationary three-dimensional MHD equations with held-aligned incompressible flows 
are calculated. Other applications of the transformation properties as a useful method for 
hnding explicit solutions for the general MHD equilibrium equations starting from known 
solutions can be found in RefsJSHSl whereas in RefJSIthe equivalence is shown of this method 
with the classical approach to infinitesimal Lie point symmetries. 

The symmetries of interest in this analysis are described by continuous Lie groups of 
transformations and we refer e.g., to the books [3-[9| for all details about this method and 
for some of its typical applications in different physical problems. It can be remarked that 
the presence of a symmetry has twofold applications: i) it makes it possible to obtain new 
solutions starting from a known one, ii) it makes it possible to look for solutions which are 
left invariant by the symmetry transformation (this may be obtained considering suitable 
reduced equations). 

The analysis of the symmetry properties of the specific GS equation that describes the case 
of stationary, azimuthally symmetric MHD plasma equilibria without flows (static equilibria) 
was developed in previous papers, see e.g. Refs JTnl - IT^ and Refill for the Euler equations 
with swirl. In RefJTTl an analysis of the symmetry properties of a class of partial differential 
equations of relevance for plasma physics, including the GS equation as a special case, has 
been performed. In particular, in RefJTOl the notion of “weaker” symmetries (conditional 
symmetries and similar) was introduced for the GS equation; the presence of such symmetries 
made it possible to identify for the first time some additional classes of solutions to the GS 
equation that describe D-shaped toroidal plasma equilibria. The present paper is devoted 
to the investigation of a generalized form of the GS equation that describes two-dimensional 
azimuthally symmetric plasma equilibria in the presence of poloidal and toroidal plasma 
flows, see Refs JTSl - ITTl but we will restrict our study to incompressible flows. The rationale 
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behind the incompressibility assumption stems from the fact that it has been known in the 
literature, RefsJTHH22l that in this limit the equation that takes the role of the GS in the 
presence of plasma flows, Eq.Q in Sec|n| can be transformed into a much simpler equation 
of the GS-type. This is made possible by a redehnition of the dependent variable (the 
magnetic flux function -0) and can be applied as long as the poloidal flow is sub-alfveni(p3!. 
This method allows us to generalize to a plasma with flows a class of solutions identihed in 
RefJTOl using this symmetry approach. 


The present paper is organized as follows. In SecjlTjwe introduce the Generalized Grad 
Shafranov (GGS) equation which includes the effect of poloidal and toroidal flows, take 
the incompressible plasma limit, dehne the new dependent variable y = x{'4’) c^nd obtain a 
modihed GS equation with an additional term that has an dependence of the cylindrical 
radius r. 

In Secjllll we apply the concept of weak conditional symmetry to the GS equation. First we 
recover the D-shaped solution of the GS equation described in RefJTOl as a special case of 
a wider class of solutions. The common geometrical structure in the poloidal r-z plane of 
this class of solutions is discussed, and shown to arise directly from its symmetry properties. 
On the other hand, determining the specihc spatial shape of these solutions involves the 
(numerical) integration of a nonlinear second order ODE whose coefficients depend of the 
physical parameters of the specific equilibrium configuration. A qualitative classification of 
the different types of solutions of this ODE is given. 

In Sec TV A we extend this analysis to the GGS equation and describe a novel class of equi¬ 
libria solutions with flows. These solutions share the same common geometrical structure of 
the solutions without flows because the symmetry properties are the same. Specific exam¬ 
ples of equilibria with flows are then shown, including solutions that exhibit three magnetic 
axes along the toroidal direction. These latter solutions correspond to configurations made 
of three separate sets of tori, each set nested around one of the magnetic axes. By imposing 
appropriate boundary conditions corresponding to '0='0=consf surfaces, each of these sets 
can be considered as an independent equilibrium configuration. 

In the Gonclusions we note that in the case where poloidal flows are present the reconstruc¬ 
tion from the solutions of the GGS equation of the explicit spatial dependence of the plasma 
pressure, of the poloidal and toroidal velocity fields and of the toroidal field from that of the 
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variable x involves the inversion of the integral relationship between ^|J and x together with 
the solution of a set of nonlinear algebraic equations. This inversion is not required in the 
case of toroidal flows and an explicit example of the effect of the toroidal flows and their 
gradients on the equilibrium conhguration is shown. 

Finally, we indicate further extensions of the solution procedure of the GS and of the GGS 
equations employed in the present article. 


II. GENERALIZED GRAD SHAFRANOV EQUATION 


We consider an axisymmetric Magnetohydrodynamic plasma equilibrium with poloidal 
{vp) and azimuthal (v^) flows and write the corresponding generalized Grad Shafranov equa¬ 
tion (see Refs JT31 - [T7|) in the form 


V- 


1 - 


47rp 




+ 


47rp 


= -AirpiJ' -UsT’ + rv^g') 


-^{n + rv^T)in' + rv^r). (i) 

Here the notation of RefJTTI has been adopted with coordinates (r, (p, z) where the azimuthal 
angle cp is an ignorable coordinate. In Eq.([^ a prime denotes differentiation with respect 
to the magnetic flux function g, p{r, z) is the plasma density and the magnetic held B is 
expressed in terms of its azimuthal and poloidal components, and Bp through the hux 
function '0(r, z) as 


B = + Vg X Vp , 

(2) 

where p = rVp is the unit vector in the azimuthal direction. 

are dehned as 

The hux functions in Eq.([^ 

7 (V’) = Anpvp/Bp, 

(3) 

ggp) =v^/r - 7Bp/{A7rpr), 

(4) 

Ti (g) = rBp - r7v^, 

(5) 

iig) = s, 

(6) 

J (pP) = u + pUp + v^/2- rVpQ. 

(7) 


where S is the plasma entropy, U is the plasma internal energy ad Up its derivative at con¬ 
stant entropy. 
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In the following we consider an incompressible equilibrinm with p = const, in which case 
p{J'' — UsT') = {p + pv‘^/2 — prv^Q)' and dehne the Alfvenic Mach nnmber M{il>) = 
J-/{4:TipY^^. Inverting Eqs.(|4jj^ we obtain 

rB^ = [r‘^Tg + n]/{l-MY, ( 8 ) 

rv^ = [r^g + W/(47rp)]/(l - M^). (9) 

The r.h.s. of Eq.Q contains fnnctions of Y and powers of r. By substituting for rv^p its 
expression from Eq.([^ and collecting powers of r^, we can rewrite the r.h.s. of Eq.([^ as 


1 




1-M2 


gj^n 


+ A7Tp{J' -UsT) 


—r 


{g^j^V2y 

1 - M2 


+ 47rp(^V2)' 


( 10 ) 


1-M2^ 

In the limit of zero plasma flow the 1/r^ term reduces to the toroidal magnetic held 
contribution to the standard GS equation and the r-independent term to the plasma pressure 
contribution. The plasma poloidal how modihes the operator on the l.h.s. of Eq.([^ as 
discussed in RefsJT^ and [161 On the r.h.s. the poloidal how introduces the multiplication 
factors depending on (1 — M^)“^ and both hows modify the hux functions in the 1/r^ in 
the r-independent term. In, addition they introduce a new term proportional to r^. In 
the remaining part of this article we will consider only sub-alfvenic poloidal hows so that 

m2 < 1. 


A. A new dependent variable: y = xiY) 


The fact that the alfvenic Mach number M in a constant density equilibrium is a hux 
function allows us to bring back the operator on the l.h.s. of Eq.([^ to its standard GS form 
by redehning the dependent variable Y- Following Refs JT8l - l^ we dehne the new dependent 
variable y = y[t/’(r, z)] as 


x{Y)= / [1 - m 2 ( p )] 1 / 2 ^ p , 


( 11 ) 


where < 1. Then, we obtain 


^2 


dvrp J r2 




dvrp 7 


= (1 - M2)l/2y. 




( 12 ) 


Finally, re-expressing the hux functions of Y on the r.h.s. of Eq.(lO) as functions of y 


through their dependence on ijj{x) ^md using cylindrical coordinates, we obtain (see also 
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Refj2l] and |25]) 


(13) 


where 


Idx d^x A ^ \ < 2 y, / N , i A ( \ 
ITT “ + ITT = -^oix) + r A2 {x) + r A(x), 

or^ r or a 


A(x) = 


A(x) = 


A(x) = 

1 

[1 - M2] 1/2 
1 


2[1 -M2]i/2’ 


V? 


1 - M2 


Q7U 

1 - M2 


[1 -M2]l/2 [ 1 - M2 


+ 47rp {J' - UsX') 
+ 47rp(6;V2)' 


(14) 

(15) 

(16) 


III. SYMMETRIES, CONDITIONAL AND “WEAK” CONDITIONAL 
SYMMETRIES. 


A. GS equation with no flow: Aa = 0 


Let us start by considering the “standard” GS equation (i.e. with no flow so that in 
particular x('0) = ■?//) 

rfli!) 1 rill) rfli!) 

(17) 




Qj.2 ^ Q^2 

The classiflcation presented in Ref JTnl exhausts all the possible Lie point-symmetries admit¬ 
ted by this equation. The presence of these symmetries is strictly related to precise choices 
of the functions F{'ip) and as fully discussed in RefJTOl Other papers devoted to the 

symmetry analysis of this equation (or of particular cases of it) are e.g. RefsJTTHTH In 
Ref|2Sl also an example of generalized symmetry is proposed. 

As shown in RefJTUl one can also look for weaker notions of symmetries, and in particular 
for the existence of conditional symmetries Refs l271 - I2ni (see also RefsHH] and EDI for other 
references). Let us recall that a conditional symmetry Y generates a transformation which 
does not map solutions into solutions, but defines a Y-invariant variable with the property 
that the solutions of the reduced equation obtained writing the initial one in terms of this 
variable, are also solutions of the initial equation. 

In Ref. [ini, we have shown to be useful to introduce even weaker type of conditional 
symmetry (see Refjrni [2ni and EU] for this notion of symmetry); in particular, the choice 


d d 

or Oz 


(18) 
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has revealed to be especially convenient (see RefJTOl for a discussion and other details). 
Writing the GS equation in terms of the F-invariant variable 


s 


— 


one obtains the following equation, where the variable r still appears (playing here the role 
of a parameter, as happens in the case of weak conditional symmetries) 

(8r^ - 4:s)^jJss - + Gii/j). (19) 

This equation can be naturally split into two ODE’s involving only the new variable s: 


S'll^ss = F{%Ij) and dsi/'ss + 2?/^^ = - ^(-0) • (20) 

Clearly, these equations admit solutions only if there is a precise relationship between E('0) 
and G{'ip). A simultaneous solution of these equations can be found of the form 

^{s) = (r^ - ^2)--? (21) 


for any real q ^ 0, with 


F{^jJ) = 8q{q + 


G{i/j) = -2g(2g + 1) • 


( 22 ) 


Let us notice that with the choice (22) for F, G, the GS equation admits the Lie point 
scaling symmetry 

Xi = r-^ + z-^-2q'il)-^ (23) 

or oz otp 


and that the above solution (21) is invariant under this symmetry. 


B. A D-shaped equilibrium solution. 


The choice q = —1/4 in (21) is of special interest as it makes it possible to construct a 


new class of solutions using an exceptional symmetry of the GS equation. We consider the 
solutioii^ 

^ ^ fiiA (24) 

which clearly holds only in the region \z\ < r, and solves the GS equation with 


F{^l^) = -(3/2)^-^ , G(^) = (1/4)^ 


-3 


(25) 
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As shown in RefJTOl with the above choice for F, G and with q = —1/4, the GS equation 
admits also the “exceptional” symmetry (see also Ref ill 

X, = + + ( 26 ) 

The presence of this symmetry implies that if ■^(r, z) is a solution of the GS equation, then 
also 

^ir,z) = [C{r, z, X)Y^^^p{r{r, z, X),z{r, z, X)) (27) 


where 

C{r, z,X) = 1 + X\r^ + z"^) + 2Xz 
r = r[C{r,z,X)]~^ and T= [z + X{r‘^ + z'^)][C{r, z, X)]~^ 
is a solution of the equation for any value of the real parameter A. Taking advantage from 


this symmetry, we can construct from the solution (24) a continuous family of solutions (see 
Refiini): we obtain, for any A, 


'ilj{r,z) = [r^ — (^z + X{r^ + z^)yY^^ (28) 

which holds in the interior of the two circles centered resp. in ro= 1/(2A), Zo= — 1/(2A) and 
To = —1/(2A), 2 ;o = ~1/(2A), both of radius l/(-\/2A), excluding their intersection (it is not 
restrictive to assume A > 0). Thanks to the obvious invariance of the GS equation under 
translations z ^ z + const., we can shift this solution and choose zq = 0. The resulting 
“D-shaped” solutions and conhgurations are shown and fully discussed in RefJTOl 


C. Additional “exceptional symmetry” solutions. 

It is possible to generalize the D-shaped solution described above by taking 

F(V’) = , Gi'ip) = ao'4)~^ (29) 

where Oq and 02 are free real coefficients and by setting 

il){r,z) = y = 'r‘^/s>l. (30) 


The variable y is an invariant variable under the symmetry ( [2^ and is a trigonometric 
function of the “latitude” angle 6, (1 — 1/y = z^jr'^ = tan^ 6^) with y = 1 corresponding to 



the equatorial (z = 0) plane. For the sake of dehniteness we restrict ourselves to positive 
values of Oq and conveniently rescale the variable "0 such that for all cases Oq = 1/4. Inserting 
Eq. (30) into Eq.([l7l) one obtains, thanks to the invariance of Eq.([I^ under the scaling 
symmetry Xi ( [2^ , an O.D.E., which is given by 

(8y^ - 12y^ + 4y) 0" + (12y^ - lO?/) 0' + (-3y/2 + 1/4) 0 = (31) 

= 1 /( 40 ^) + a2y/(t)^. 


For 02 = —3/2, Eq.(31) is solved by 0 = 1 and returns the result of Sec, III B For different 


values of 02 , keeping 0(1) = 1 as a normalization condition on 0, Eq.(31) can be easily 


solved numerically. Note that the exceptional symmetry (26) depends only on the scaling 


of E(0) and of (7(0) on 0 and not on the specihc values of the coefficients oq and 02 (see 
also Sec. IVA). Then the transformation (27) applies also to the solutions of Eqs.( pOpT ) 
and sends solutions into new solutions of the form 


0(r, z) = [r^ ~ (^ + 0(r^/[r^ ~ + A(r^ + ^^))^]), 

which can be rewritten in a geometrically more transparent way as 


(32) 


(33) 


where we have made explicit the dependence of the solution on the parameter 02 , and 

yx = r'^/[r'^ — [z + X{r‘^ + z'^)Y], 1/|/a = 1 — (tan6* + Ai?/cos6*)^, (34) 

with + z"^. Notice that y = const, corresponds in the r > 0 half-plane to the 

lines r = k\z\, with k > 1, whereas yx = const, corresponds to the borders of the lunes 
(crescents) dehned by the circles 


(r q: (1 - 1 /j/a)‘''V(2A) = (2 - 1/!/a)/( 2A)^ and r > 0, 


(35) 


as shown in Figj^ (here the variable 2 ; has been translated, see below Eq.(28)). Note 
that the shape of these lunes does not depend on the value of 02 and that the function 
(l>{y\iO> 2 )/y]!^ is constant on their borders. Thus the shape in the r-z plane of the class 
of solutions discussed below can be visualized as given by the level curves of times a 
function that depends on 02 and that is constant on the borders of the lunes in Figjl] whose 
shape depends only on the transformation parameter A. 
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Figure 1: Lunes with borders given by the circles in Eq.(35) for A = 1 and 1 — l/yx = 1, red, 
1 — 1/?/a = 1/2, green and 1 — I/t/a = 0, black. The red border is the transformed of the lines 
\z\ = r, the black one of z = 0. The two (common) vertices of the lunes are the transformed of 

the origin and of the point at oo. 


Noting that the coefficient of (p” in Eq.(31) vanishes sX y = 1 and restricting to solutions 
where cp" does not diverge, keeping as mentioned above 0(1) = 1, we hnd that the value of 
0' at ?/ = 1 is given by 

2 0'(1) =3/2 + 02. (36) 


Two classes of solutions of Eq.(31) can be identified depending on the value of 02 - For 
02 > —3/2, and then 0'(1) > 0, the function 0 remains positive over the whole interval 
1 < y < 00 and behaves asymptotically as <p{y) oc y^^"^ for y —)■ cx) i.e., for 2 ;^ —)■ . Note 


that for 02 = —1, one has the exact solution <p = y^^'^ of Eq.(31), which yields the special 
solution 0(r, z) = 

For 02 < —3/2, the function <p{y) reaches zero at a finite value y = 1 /( 02 ). For y ^ y, 
i.e., for z'^ —)■ r^(l — l/y), it behaves locally as (p{y) oc {y — yY^"^ which corresponds to 
YiT^z) OC — r^(l — l/yPY^^. In this case a complementary solution can be found for 
y < y < 00 which starts as (p{y) oc {y — yY^^ and behaves as oc y^^^ for ?/ —)■ 00 . 

The structure in the r-z half-plane of the solutions with 02 > —3/2, after being transformed 


according to Eq.(27), is similar to that of the D-shaped solution obtained in RefJTnl with the 
magnetic axis shifting towards larger values of r with increasing values of 02 - The structure 
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of the solutions with 02 < —3/2 will be discussed in the next section in the frame of the 
solutions of the GGS equation. 


IV. EQUILIBRIA WITH FLOWS 
A. Symmetries of the GGS equation 


We now consider the case of the generalized GS equation (13). 


It is easy to show that the two symmetries Xi and X 2 given above (23) and (26) are still 


valid: more precisely, the symmetry Xi is admitted by any equation of the form 


d^X _ ^ 

Qj.2 j. Qj. Q^2 


E l 1 _|_ 2 +^ 

air X "" 


(37) 


and, similarly, the exceptional symmetry X 2 is admitted by any equation of the form 


d'^X _ ^ 

dr 2 r dr dz^ 




(38) 


where at the r.h.s. any combination of terms r^x^^ (resp. r^x ^ with arbitrary 
constants ai and - in principle - any real value of ^ is admitted 


Restricting to the cases relevant to Eq.(13), we can write instead of (37) 


d‘^X , d'^X 

dr‘2 


r dr dz"^ 


1+1 


2 . i+; 




= aox 5 + asr X -7 + a 4 r y 9 


and hnally, with q = —1/4, 


d'^X 1 dx d‘^x ®o 02+^ 

1“ v: v Q "h ^ h 7 
X X X 


d '\/H 


(39) 


(40) 


B. Solutions of the GGS equation 


We proceed as in Sec. Ill G by dehning 


X(r, z) = 0 ( 1 /), y = r^/s > 1, 0(1) = 1, 


(41) 


and obtain an O.D.E. (thanks again to the invariance of Eq.(40) under the symmetry X 


which is of the form of Eq. (31) with an additional nonlinear term (as before we put oq = 1/4): 


(8/ - 12/ + iy) 4 ," + (12/ - lOj/) / + {-3y/2 + 1/4).(. = 


(42) 
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= 1/(40-^) + 02 y/cjy + 04 


As in Sec. (III C) the exceptional symmetry allows us to construct a family of D-shaped 


equilibria!^ which now include plasma flows 


(43) 


These solutions can again be visualized as the level curves of times a function that 
now depends on both 02 and 04 and that is constant on the borders of the lunes dehned in 
Sec JIITCl 


Note that for 04 7 ^ 0 there is no constant solution 0 = 1, while Eq.(36) is changed into 

2 0'(1) = 3/2 -|- 02 + 04. 


(44) 


Three classes of solutions of Eq.(42) can be identihed depending on the values of 02 and 
04 as indicated in Figj^ 

•* 04 



02 


Figure 2: Regions defining the three different types of solutions of Eq.(42) in the 02 and 04 


parameter plane. For 04 = 0 see also Sec III C 


For 0'(1) = (04 -|- 02 + 3/2)/2 > 0 with the additional condition indicated by numerical 
calculations 04 -|- p(o 2 + 3/2) > 0 with p ~ 0.85, region A of Figj^ the function 0 remains 
positive monotonically growing over the whole interval 1 < y < 00 and behaves asymptoti¬ 
cally as </)(y) oc y^!^ for ?/ — )■ 00 as in the case without flow. 

For 04 -|- p(o 2 -|- 3/2) < 0 and 04 < 0, region B of Figj^ the solution reaches zero at a hnite 
value y = y, as shown in Figj^in terms of the variable 1 — 1 /y = In this case, 0(?/) 

behaves locally as 0 (|/) oc {y — yY^^ and, as in the case without flow, a complementary 
solution can be found for y < y < 00 which starts as 0 (j/) oc [y — yY^^ and behaves for 
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1.4 


(/.(y)/i/V4 



1-1/y 


Figure 3: Plot of (/>(y, 0.5,-1.85)/y^/^, see Eq. ( |43[ ) , expressed as a function of the variable 1 — 1/y 
whose value determines the lune borders in Figj^ Part of the solution for y > y is not shown. 





i-i/y 

0 

Figure 4: Plot of cj){y, —3.5,0.5)/y^/^, versus the variable 1 — 1 /y. 


?/ —)■ oo as oc 

For 04 + 02 + 3/2 < 0 and 04 > 0, region C of Figj^ we have a solution which remains strictly 
positive but does not behave monotonically, and asymptotically grows as oc as shown 
in Figi in terms of the variable 1 — 1/y. For 04 + p(o 2 + 3/2) ~ 0, there are “separating so¬ 
lutions” which behave asymptotically as Notice the particular case 04 = 0, 02 = —3/2 

where one recovers the known solution 0 = 1 . 

In order to illustrate the different shape in the r-z plane of the solutions belonging to 


region A and to regions B-C, when transformed according to Eq.(43) and shifted in z, we 


may use simple analytical htting formulae of the numerical solutions of 0 (?/, 02 , 04 ). 
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1.4 



Figure 5: Shaded contour plot of the solution - 2 ) = “l-Sj for A = 1. 


The shaded contour plots of two such solutions are shown in Figsj^ and over to whole 
domain within the red lune in Figjlj Obviously, these solutions can be restricted to within 
a reduced domain inside this lune by imposing (conducting) boundary conditions on the 
y = const. (—)■ V’ = const.) border of the chosen domain. 

The solution from region A with 02 = —1.8, 04 = 0.7 and A = 1 is shown in Figj^ it is 
D-shaped with its magnetic axis shifted toward large r and steep gradients of y near the 
border of its outermost magnetic surface. 

The solution from region C with 02 = —3.5, 04 = 0.5 (the same of Figj^ and A = 1 is shown 
in Figj^ It exhibits three magnetic axes along the toroidal direction. Thus this solution 
corresponds to a conhguration made of three sets of surfaces, each set nested around one of 
the magnetic axes. Each set can be considered as an independent equilibrium configuration. 
The configuration centered at the middle axis is D shaped while the outer configuration has 
a crescent shape. A similar structure, but with steeper gradients, is obtained for solutions 
from region B. In this latter case the three sets of surfaces are separated by the y = 0 
surface. 
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Figure 6: Shaded contour plot of the solution - 2 ) = “3.5, 0.5)/?/]^^^ for A = 1. 


V. CONCLUSIONS 


In this article we have developed a method based on the Lie symmetries of the GS and 
of the GGS equations that we have identified by imposing specific relationships between the 
flux functions of the magnetic flux variable in these equations. 


In the case of the GS equation (plasma equilibria without flows) we have extended the 
class of the D-shaped solutions found in RefJTUl Given the numerical solution of an ODE 


whose general properties have been discussed in Sec jlll G[ these solutions can be worked out 
explicitly in order to obtain the dependence of the plasma pressure and of the toroidal field 
on the cylindrical coordinate r and z as done in RefJTOl 


In the case of the GGS equation, when treating equilibria with poloidal flows, we have 
made use of an implicit integral transformation from to a new dependent variable x{'4^) 
that can be applied when the alfvenic Mach number of the poloidal flow is a flux function 
and is smaller than unity. Although the shape of y (and thus of V’) can be found with a 
relatively minor generalization of the procedure adopted for the GS equation, finding the 
explicit dependence of the plasma pressure, of the poloidal and toroidal velocity fields and 
of the toroidal field on r and z in the GGS case requires the solution of a set of nonlinear 
equations and, in the case of poloidal flows, the inversion of the integral relationship between 
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0 and X once the flnx fnnction -7^(0) is chosen. The system of Eqs. ( p!4p!5pT| ) must in fact 
be solved for the remaining flux functions to make the functions ^ 2 , *^4 compatible with 


Eq.(39). 



- 1.5 -1 - 0.5 0 0.5 


z 

Figure 7: Shaded contourplots of 0(r, z) for A = 1, 02 = —3/2 and 04 = 0 top frame, 04 = 0.1 

center frame and 04 = — 0.1 bottom frame. 


This aspect of the problem, together with possible restrictions it might bring on the physical 
range over which the coefficients 02 and 04 can vary, has not been addressed explicitly in 
the present article. However if physical boundary conditions are imposed, for example by 
external conductors, that restrict the solution to the “flat” central domain of e.g., the con- 
hguration shown in Figj^ that is to say to a domain inside y = const (and thus 0 = const) 
surfaces within which y is sufficiently larger than zero and the Mach number sufficiently 
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smaller than one, Eq.(p!T| can be inverted by a simple pertnrbation procednre. On the con¬ 
trary, the inverse power dependence on y of the terms of the terms on the r.h.s. of GGS 


eqnation (40) implies that the eqnilibrium qnantities, including the flow prohles, will develop 
large gradients, and/or the Mach number will approach unity, at the borders of the lunes 
in Figj^ where x 0 (and, for solutions in region B, near the internal x = 0 curves). This 
can be seen by combining d/dilj = [1 — d/dx with the x derivatives of the flux 


functions ^o(x))“^ 2 (x) * 4,4 (x) in Eqs.(14 16). 

Furthermore, it is instructive to compare the D shaped equilibrium conhguration without 
flows discussed in Ref JTnl with two corresponding equilibria (at hxed J'i'ip), see Eq.([^) with 
toroidal flows but without poloidal flows so that the functions x and ip coincide. This 
comparison is shown in Figj^ wher^^ the shaded contour plot of the D-shaped equilibrium 
of RefJTUlis shown in the top frame (02 = —3/2, 04 = 0 ), that of a corresponding equilibrium 
with angular velocity increasing with ip (02 = —3/2, 04 > 0) is shown in the center frame 
and one with angular velocity decreasing with ip (02 = —3/2, 04 < 0) in the bottom frame. 
We see that a positive gradient of the centrifugal term due to the toroidal flow (04 > 0) 
leads to a shift of the maximum of the flux function ip{r, z) towards larger values of r and 
to the formation of large gradients at the outer border, as characteristic of the equilibrium 
solutions in region A (see e.g., Figj^. On the contrary, a negative gradient of the centrifugal 
term (04 < 0 ) causes ip{r, z) to vanish not only at the boundary of the solution domain but 
also along two curves inside the domain as is the case for the solutions in region G. This 
splits the solution into “separate” solutions with different magnetic axes, as characteristic 
of the equilibrium solutions in regions B and G (see e.g., Figj^. From the bottom frame 
of Figj^we see that the inner solution is, for the chosen parameters, vanishingly thin, the 
central one has a rather flat radial prohle while the outer solution has large gradients. 

As a concluding remark we observe that the procedure developed in the present paper 
can be further extended to a wider class of solutions by adapting it to the different choices 
of the variable s = s(r, z) made in Sec.IV of RefJTOl 
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